syms s;
P=(s+3)^2*(s^2+3*s+2)*(s^3+12*s^2+48*s+64)
P1=simplify(P),
P2=expand(P),
P3=factor(P),
P4=prod(P3)